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bJ)' Abstract 
^ . 

. Differential-difference integrable exponential type systems are studied corresponding 

to the Cartan matrices of semi-simple or affine Lie algebras. For the systems corresponding 

to the algebras A2, B2, C2, G2 the complete sets of integrals in both directions are found. 

For the simple Lie algebras of the classical series Ai\f, Bjsf, Cm and affine algebras of series 
(2) 

. the corresponding systems are supplied with the Lax representation. 

■ Keywords: semi-simple Lie algebra, affine Lie algebra, differential-difference systems, x- 

integral, n-integral, characteristic Lie algebra, Toda field theory. 

(M ■ PACS number: 02.30.Ik 

^ 1 Introduction 

10 '. It is well known (see [1]) that one can assign an integrable exponential type system of hyperbolic 
^ ' equations to each semi-simple or affine Lie algebra. Denote through A = {ajj} the Cartan 
matrix, corresponding to the Lie algebra, then the desired exponential system is written as 
follows 

rly = e^'''^'' . (1) 

I The classical Toda field theories in the continuous space-time are based on infinite dimensional 
integrable Hamiltonian systems closely connected with ([1]). In this context these equations have 
been intensively studied in mathematical physics during the last three decades p[^-|12j. 

In the present article we study the problem of discretization of integrable systems of the form 
([1]) preserving integrability, it is important since discretization makes a bridge between quantum 
and classical field theories. More exactly, we look for integrable differential-difference systems, 
with a discrete space-like variable n and a continuous time x which correspond to ([1]) in the 
continuum limit. The problem of discretization of exactly solvable models is intensively studied. 
For the particular case x = y Suris [T3] exhibited discrete integrable versions of the Toda chain 
corresponding to all the non-exceptional affine Lie algebras. Today different approaches to the 
problem are known in literature |13]-|20]. The approach using the Hamiltonian formulation of 
the Toda field equations is studied in |T7j, where the equation of motion for the Toda lattice. 
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associated with the simple Lie algebras, is derived. In order to approve integrability of the 
lattices suggested the authors studied the properties of the corresponding r-functions. Further 
discretization of the lattice model leads to difference-difference equations on discrete space- 
time, having a very large variety of applications in Physics, related "... to commuting transfer 
matrices in solvable lattice models, q-characters of Kirillov-Reshetikhin modules of quantum 
affine algebras, cluster algebras with coefficients, periodicity conjectures of Zamolodchikov and 
others, dilogarithm identities in conformal field theory, difference analog of L-operators in KP 
hierarchy. Stokes phenomena in ID Schrodinger problem, AdS/CFT correspondence, Toda field 
equations on discrete space-time, Laplace sequence in discrete geometry, Fermionic character 
formulas and combinatorial completeness of Bethe ansatz, Q-system and ideal gas with exclusion 
statistics, analytic and thermodynamic Bethe ansatz, quantum transfer matrix method and so 
fortlj^'. 

Our approach to the problem of constructing discrete analogues of system ([1]) is based on 
the technique of integrable cutting off constraints [19], [21] and on the integrals preserving 
method of discretization, suggested recently in [20] , for the Darboux integrable hyperbolic type 
models. 

Define a mapping assigning to each N x N matrix A a system of differential-difference 
equations of the form 

rl,.-rl = (^<^'^^^^''^^''\ z,j = l,2,...,iV, (2) 

where the given matrix A is decomposed into a sum of two triangular matrices: A = j4_|_ + A_ 
with A+ = {afj} being upper triangular and A_ = {a~} - lower triangular matrices such that 
all diagonal entries of A+ equal unity. Here the functions = r^{n,x), j = 1,...,N are the 
searched field variables. Everywhere below the subindex denotes a shift of the discrete variable 
n or the derivative with respect to x: rl = r^{n + k,x) and r-^^j = -^r^{n,x). Sometimes 

we use also = -^r^{n,x), r^^ = ^r^{n,x), ... . We refer to ([2]) as a differential- difference 
exponential type system. 

The further investigations approve the following hypothesis. 
Conjecture 1. a) If A is the Cartan matrix of a semi-simple Lie algebra then is Darboux 
integrable. More detailed explanation of Darboux integrability is given below in Definition fTH 
Roughly speaking Darboux integrable equations can be solved in quadratures, 
b) If A is the Cartan matrix of an affine Lie algebra then ([2]) is S'-integrable (it admits large 
classes of exact multi-soliton solutions). 

Part a) of the conjecture is checked for the following cases: A2, B2, C2, G2. In these cases 
the exponential system ([2]) reads as follows 

rU-rl = e^'^^'i-^l (3) 

^2 _ 2 _ -cr^+r'^+rf 
' l,x ' X 1 

where c = 1,2,3. Integrals in both directions for these systems are found. For the systems 
related to the algebras of the series Ajy, with N > 2 the Lax pairs are given. They are also 
Darboux integrable. Note that the systems corresponding to the Cartan matrices Cjv and Sat, 
N = 2, 3 ... , are the same up to renaming the variables. The system corresponding to Ajy is 
easily reduced to that found earlier in |l3j . 

*A. Kuniba, T.Nakanishi, J.Suzuki, T-systems and Y -systems in integrable systems, J. Phys. A: Math. 
Theor. 44 (2011) 103001 (146pp) 
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Part b) is illustrated by an example of system corresponding to the Cartan matrix of 

(2) 

algebra D]^ for which the Lax pair is found by using the integrable cutting off constraint for 
infinite chain of the Backlund transformations for the classical two-dimensional Toda lattice. 

Let us recall the necessary definitions (see [6]). Function w = t(;(r, r[i], r^j, rjfc]) with 
r = (r^, r^, ...r^), r[j] = |^ is called a y- integral of order k of system ([1]) if its total derivative 
with respect to y vanishes DyW = 0. Similarly, function iv = ?I;(r, f[i], f[2], f[fc]), where 
f[j] = 1^ is an x-integral of order k if D^w = 0. System ([1]) is called Darboux integrable if 
it admits ?y-integrals w^^\ w'^'^\ . . . , w^^^ of orders q^^\ q^'^\ . . . , g*^^-* correspondingly and N 
x-integrals w^^\ 'w'^'^\ . . . , w^'^^ of orders . . . jp*-^-* such that 

det\——]^Q and det U 0, A:, j = 1, 2, ...TV. (4) 

In this case the set of 2N above integrals constitute a basis of integrals of the system ([T]). Inte- 
grals w'^^\ w^'^\ . . . , w'^^'> (or integrals w'^^\ w^'^\ . . . , w^'^^) satisfying are called independent. 

Analogously we define integrals for the differential-difference chain ([2]). Here the set of 
dynamical variables consists of the field variables r^, r^, . . . , r^, their shifts D'^r^ = rj^ and 
derivatives -D^r-' = r^.j. Recall that D is the shift operator acting as follows Dh{n) = h{n + 1) 
and Dx is the operator of the total derivative with respect to x. For the convenience we put 
r = (r^, r^, . . . , r^). Function F = F{x, r^', rk'+i, ■ ■ ■ , r^), where k' < k is called an x-integral 
of system ([2]) of order k — k'ii D^F = 0. And / = I{x, r, r[i], r^j, . . . , r^^]) is called an n- integral 
of order k if DI = L 

Definition 1. We call system ^ Darboux integrable if it admits N x-integrals F^^\ F^'^\ . . . , F^'^^ 
of the form F^''^ = F^^\x, r, ri, . . . , r^) of order k and N n-integrals . . . , J^^-* of order 

s such that 



det[-—j^O and det \^-^ ] ^ 0, ^, j = 1, 2, . . . A^. (5) 



A set of 2A^ integrals satisfying ([5]) constitutes a basis of integrals of system ([2]). Integrals 
pW ^ ^(2)^ _ _ _ ^ p(N) ^^Yl as integrals . . . , J*-^-' satisfying ([5]) are called independent 

because non of them can be functionally expressed through the other integrals and their shifts 
and derivatives. 

Conjecture 2. Suppose that the function 

w = w{r, rx,rxx, ...) 

is a ^/-integral for the system ([T]) with the Cartan matrix ^4 of a semi-simple Lie algebra. Then 
this function is also an ri-integral for the system (|2]) with the same matrix A. 

Conjecture 2 is approved by the examples below (see section 3). For instance, the functions 



= w^'^=rl + rl-irir + rlrl 
= ^(2) ^ ^1 _^ ^ic^2 _ ^ _^ c^l^2^2 _ ^1/ : 



XX ' XX \ Xj 'XX \ xl 1 

1^»,2\2 



are y— integrals for the system 



1 2r^-r^ 

r = p 

' xy ^ ' 

^2 ^ -rl+2r2 
' xy ^ 
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and n— integrals for the differential-difference chain 



^1 -r^ = f,r^+rl-rl 
' l,x ' X ^ 1 

rl-rl = e-^+"^+^?. (6) 

A set of independent x-integrals for chain (Q is given in Example [L] 

To show the existence of x- and n-integrals one can consider characteristic x-algebra and 
characteristic n-algebra L„ (see [22] and references therein). 

Theorem 1. (see ) A system of equations admits a non trivial x-integral if and only if its 
characteristic x-algebra is of finite dimension. 

Theorem 2. (see 122]) A system of equations admits a non trivial n-integral if and only if its 
characteristic n-algebra is of finite dimension. 

Conjecture 3. Suppose that the evolutionary type system 

fiX^ ^Xy "^XXI ■ ■ ■ ) 

is a symmetry of system ([1]) with the Cartan matrix of a semi-simple or affine Lie algebra. 
Then this equation is also a symmetry of system ([2]) with the same matrix A. 

The paper is organized as follows. In section 2 we introduce characteristic x-algebra for a 
system of the form ([7]) (see also ([3])) and prove that for the values of c = 1, 2, 3 the algebra is of 
finite dimension. This allows one to deduce a system of linear first order PDE the x-integrals 
should satisfy to. These systems of equations are successfully solved and independent integrals 
are found for each of the cases c = 1,2,3. In the third section n-integrals for system ([?]) are 
presented. The close connection of the integrals found with ^/-integrals for the corresponding 
continuous systems (II]) is demonstrated. 

In fourth section an alternative approach is considered to construct finite field systems. In 
this approach one uses the cutting off boundary conditions preserving integrability for the open 
chains (see |19], |2TJ). By using certain special choices of cutting off constraint we have been 
able to reproduce in this way some of the chains obtained by applying the mapping above (see 
(E])). The second method to construct integrable systems is also important since it allows one 
to find the Lax pairs for the finite field reductions. 



2 Characteristic x-algebras and x-integrals 

Let us consider the problem of finding the x-integrals for the system ([3]) which is rewritten for 
the convenience in the following form 

ai,x-a, = e'^+'^^-^\ (7) 

where c = 1,2,3, a = r^, b = r^. First we introduce the notion of the characteristic x-algebra 
for the chain ([7]). According to the definition an x-integral F should satisfy the equation 
DxF[x, a, b, ai, bi, a_i, b_i, ...) = which is equivalent to the equation ZF = where 

k=—oo 
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Another important observation is that x-integral cannot depend on the variables a^, despite 
the coefficients of the equation ZF = depend on these variables. That is why we have to put 
two additional equations YiF = 0, Y2F = where 

d d 

Ottx OOx 

Denote through the Lie algebra generated by the operators Z, Yi, Y2 over the field of 
constants. We call this algebra characteristic Lie algebra in x-direction or x-algebra (see also 
p2j). The purpose of the present section is to give the complete description of the x- algebra 
for system ([7]). 

For the system (JTj) we have 

k k 

ak,x = ax + J2 6"-^+'^'-'' and a.k,x = a.-J2 e"-+'^-'+i-''-^+i , 

i=l i=l 
k k 

bk,x = bx + Yl 6-™'-^+'''-^+^' and b_k,x = - J] 6'^''-^+^-^+^-^+' . 
1=1 1=1 

So, the vector field Z takes a form 

d d . d 

ox oa oh 

00 / k \ r) °° f ^ \ Pi 

k=l \ 1=1 / k=l \ i=l / 

00/fc \ p. 00 / k \ Pi 

+ E + E ^ + E - E e----^--'- U^. (9) 

k=l \ 1=1 J ^ k=l \ i=l / ^ 

Taking commutators of the vector fields Yi, Y2 and Z we get 

k=—oo k=—oo 

Taking commutators of the vector fields A, B and Z we get vector fields A* = [A, Z] and 
B* = [B,Z]: 

00/fc \rjOo/fc \ f) 

fc=i V i=i J k=i \ i=i J ^^-^ 



00 



+ E I E ^-^-^^ ) ^ + E h E Uf^ (10) 

fc=l V 1=1 / ^ k=l \ i=l 



and 



^* = E - E ^"-^^""'^ + E E e"--^'^--^-- ) ^+ 

fe=l \ i=l / k=l \ i=l 



00 
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It is convenient to introduce the following vector fields 



4 — c 



fc=l \i=l / k=l \i=l / 



and 



P2 = (A* + 2B* 

4 — c 



= E E ^-^^^^ ^ + E -2 E e-----+--^ (13) 

k=l \i=l / ^ k=\ \ t=l / ^ 

Obviously the vector fields 

^, Fi, F2, ^, A and P2 

belong to the x-algebra and vector field Z can be obtained as their linear combination. We use 
the following lemma to show that x-algebra is of finite dimension. 

Lemma 1. (see f2W ) Suppose that the vector field 

d 



^ = E(«^^+«-^)+E(/^^J^+/5- 

k=l ^ " '^^ k=l ^ 



(14) 



satisfies the equality DKD ^ = hK , where h is a function depending on shifts and derivatives 
of variables a and b, then K = Q. 

Let us consider the sequence of commutators Tq = P2, 71 = [Pi, P2] and T„ = [Pi,T„_i], 
n = 1,2 . . . . To apply the lemma 1 we need to know how the transformation D{*)D^^ acts on 
the vector fields A, B, Pi and P2. One can easily check that 

DAD-^ = A, DBD-^ = B. (15) 

Lemma 2. 



DPiD-^ = Pi- e'^+'^i-^'iA, 
DP2D-^ = P2- e-="+^+^iP. 



(16) 



Proof. Let DPiD'^ = yUfc- h fi-kT: , then for A; = 1, 2 . . . 



'k-l \ fc-1 



/ifc = DPiD-\ak) = DPi{ak-i) = D ( e^^-'^"^'^^ ] = ^ e"^+""+^"^^+i. 
Hence we have 



i=l I i=l 



Similar equality holds for A; = 1, 2 ... . Thus, we have DP^D^^ = Pi — e"+«i-^i^. Formula 
for DP2D~^ is proved in the same way. □ 
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Lemma 3. For any k > c + 1, where c takes the values c = 1, 2, 3 we have = 0. 
Proof. Using previous lemma we have 

DTiD-^ = D[Pi,P2]D-^ = \DPiD-\ DPaP*"^] = \Pi - e"^"^-^^ A, P2 - g-^^+^+^'ifi] = 

Ti + ce"+"^-''^P2 - 6"™+^+''^ Pi - 06^^-'^"+"'+^ B. (17) 

In the same way one can show 
DT2D-^ = T2 + 2(c - l)e''^'''-'"Ti + c{c - iya+2a^-2b^ 

- 2(C - l)e(l-^)'^+"i+^Pi - c(c - l)e^^-c)a+2a,+b-b,Q^ ^^g^) 

DT^D~^ = P3 + 3(c - 2)e°+°i-''ir2 + 3(c - l)(c - 2)e2'^+2ai-26iy^^ 

+ c(c - l)(c - 2)e3'^+3'^^-3^iP2 - 3(c - l)(c - 2)e(2-^)'^+2'^i+^-^iPi- 

- c(c - l)(c - 2)e(3-")'^+=^°i+^-2^ip, (19) 

DTiD-^ = T4 + 4(c - 3)e''+''i-''iT3 + 6(c - 2)(c - 3) g^^+s^i -2617^2+ 

+ 4(c - l)(c - 2)(c - 3)e3«+3'^i-3f'i7^_^ + c(c - l)(c - 2)(c - 3)e4'^+4«i-46ip2_ 
- 4(c - l)(c - 2)(c - 3)e(3-<=)'^+3'^i+fe-2bip^ _ c(c - l)(c - 2)(c - 3)e(4-^)'^+4ai+f>-3bi^^ (20) 

Hence the statement of the lemma directly follows from the lemma [T?] □ 

Lemma 4. For c = 1,2,3 we have [P2, Pi] = 0. 

Proof. Direct calculations show that 

D[P2,T^]D-' = [P2,T,], 

which implies that [P2,Pi] = by lemma [L] □ 

The following theorem can be easily proved by using the lemmas above. 

Theorem 3. The x-algebra corresponding to system ^ with c = 1,2,3 is of finite dimension. 

Remark 1. Emphasize that we do not study any classification problem for system ([7]) with 
arbitrary value of the parameter c, only three particular cases c = 1,2,3. The problem of 
describing all values of c for which x-algebra is of finite dimension is a subject of a separate 
investigation. 

Since the x-algebra is finite-dimensional, we can find the x-integrals corresponding to the 
values of c = 1, 2, 3 in system ([7]). 

Example 1. For c = 1 the corresponding x-algebra is generated by vector fields Yi, Y2, A, 
B, Pi, P2, Ti with the following commutators table: 





A 


Yi 


Y2 


A 


B 


Pi 


P2 


Pi 


d_ 

dx 


























Yi 


























Y2 


























A 

















2Pi 


-P2 


Pi 


B 

















-Pi 


2P2 


Pi 


Pi 











-2Pi 


Pi 





Pi 





P2 











P2 


-2P2 


-Pi 








Pi 











-Pi 


-Pi 
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The x-integral F(...a, 6, ai, 61, 02, 62 • • • ) must satisfy equations 

A{F) = 0, B{F)=0, Pi(F) = 0, P2(i^) = 0, Ti(P) 



0. 



(21) 



Solving the above system, it is enough to assume that Fi depends on b, a, bi, ai, 62, (I'^d 
F2 depends on a, hi, ai, 621 0,2, b^, we find two independent x-integrals 



P2 = e 



g-b+fei _j_ g-a+ai +61-62 _j_ gai-a2 
-a+ai _j_ gai-a2-6i+62 _|_ gb2-b3 



The system corresponds to the Cartan matrix 

2 



-1 



(22) 



and the Cartan matrix [2^) is canonically related with the simple Lie algebra A2. 

Recall that simple Lie algebra ofrankr is completely described by its Weyl generators Xj, yi, hi, 1 < 

i < r, such that 



1. Xi ^ 0, yi ^ 0, hi ^ for any i; 

2. for any i,j the relations hold 



[hi, hj 



0, 



^ij hi , 



[hi,yj] = -Nijyj. 

Here N = {Nij} is the Cartan matrix of the algebra, and 6ij is the Kroneker symbol. 

It is evident from the table of commutators that the mapping hi — )■ —A, h2 — —B, yi — )■ Pi, 
y2 — )■ P2 gives an isomorphism between the algebra generated by the operators {A, B, Pi, P2} 
and Borel subalgebra of the Lie algebra A2 generated by {hi, h2,yi,y2} ■ 

Example 2. For c = 2 the corresponding x-algebra is generated by vector fields Yi, Y2, A, 
B, Pi, P2, Ti, and T2 with the following commutators table: 





d_ 

dx 


Yi 


Y2 


A 


B 


Pi 


P2 


Ti 


T2 


d_ 

dx 





























Yi 





























Y2 





























A 

















2Pi 


-2P2 





2r2 


B 

















-Pi 


2P2 


Ti 





Pi 











-2Pi 


Pi 





Ti 


T2 





P2 











2P2 


-2P2 


-Ti 











Ti 














-Ti 


-T2 











T2 











-2T2 


















The x-integral F{...a,b,ai,bi,a2,b2 ■ ■ ■) must satisfy equations 

A{F) = 0, B{F) = 0, Pi(P) = 0, 

P2(F) = 0, Ti(P) = 0, r2(i^)=o. 
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(23) 



Solving system /[^) . it is enough to assume that Fi depends on a, &i, ai, 621 0-2) ^s? % ond 
F2 depends on b, a, hi, ai, 62, 02, 63, we find two independent x-integrals 

Fi-- 



F2 = 

In this case N 



-a+ai _j_ g-ai+a2+fe2-t'3 _j_ gai-a2-bi+62 _|_ ga2-a.3 
-b+b 



1 _j_ g-2a+2ai+6i-62 _j_ 2g^"+2ai-a2 _|_ g2ai-2a2~bi+62 _j_ g''2-63 



2 -1 
-2 2 



and we wse the map hi — )■ —5, /i2 — J- — ^4, 2/1 — P2, y2 Pi 



to identify subalgebra generated by {A, B, Pi, P2} and the Borel subalgebra of B2 generated by 
{hi,h2,yi,y2}- 

Example 3. For c = 3 the corresponding x-algebra is generated by vector fields Yi, Y2, A, 
B, Pi, P2, Ti, T2, T3 and Q = [Ti,T2] with the following commutators table: 





dx 


Yi 


Y2 


A 


B 


Pi 


P2 


Ti 


T2 


T3 


Q 


dx 



































Yi 



































Y2 



































A 

















2Pi 


-3P2 


-Ti 


T2 


3T3 





B 

















-Pi 


2P2 


Ti 





-^3 


Q 


Pi 











-2Pi 


Pi 





Ti 


T2 


T3 








P2 











3P2 


-2P2 


-Ti 

















Ti 











Ti 


-Ti 


-T2 








Q 








T2 











-T2 





-n 





-Q 











Ts 











-3T3 


T3 




















Q 














-Q 





















The x-integral F{...a, 6, ai, 61, 02, 62 • • • 
A{F) = 0, B{F) 

0, 



must satisfy equations 
0, Pi(F) = 0, P2{F) 



0, 
0. 



(24) 

Ti(F) = o, r2(F) = o, g(F) = o, r3(i^) 

Solving system [24\ ), we assume that Fi depends on a, hi, ai, 62, 02, 63, 03, 64, 04 and F2 
depends on b, a, hi, ai, 62, 02, ^3, ^3, ^4; and find two independent x-integrals 

g-a+ai _|_ g-2ai +202+62-63 _|_ 2g-'^i+2a2-a3 _|_ g-a2+a3+63-64_|_ 



Fi 
F2 



_j_g2a2-2a3-fe2+63 _|_ gai-a2-6i+62 _|_ ga3-«4 



g-6+fei _j_ 2g-2a+3ai-a2 _j_ 2g-ai+3a2-2a3 _j_ gg-a+oi +02-03 _j_ g-3a+3ai+fei-fe2 
_|_2g-a+3ai— 202-61+62 _j_ 2gai-a3-6i+ti2 _j_ g3ai— 802-261+262 _j_ g— 801+302+262- 263_j_ 
_j_3g-o+02+62-63 _j_ 2g-2ai+3o2-a3+62-63 _j_ 2g-''i+262-63 _j_ g3o2-3o3-62+63 _j_ g63- 



64 



Here the Cartan matrix is N 



The isomorphism between the subalgebra of 



the characteristic algebra generated by {A, B, Pi, P2} and the Borel subalgebra of the algebra G2 
is established by the same mapping as in the previous example. 

Remark 2. In order to find the arguments of F one has to order the set of the dynamical 
variables {oj, bi} in an alternating way as follows 

. . . a_i, b, a, bi, ai, 62, 02, 63, «3, h, 04, . . . 

and then take a segment of the length greater by one than the dimension of the system of 
equations for F. If for one of the integrals the segment begins with a then for the other one it 
begins with b. 
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3 Construction of n-integrals for differential- difference 
systems 

It turns out that due to the special form of system ([7]) it is reasonable to look for n-integrals 
using the definition. We will search an n-integrals which are in a sense homogeneous polynomials 
of the derivatives of the variables a, b. That is they are linear combinations of monomials 

"[i]"[2] • • • "[fc]"[i]"[2] • • • "[k] y^^^J 

for which the sum is constant. We also note that one can always find independent 

integrals that are linear in the highest derivative present (see p3]). Hence, to find an n- integral 
one has to fix the order of the highest derivative. Then write the general homogeneous, in the 
above sense, polynomial / with undetermined coefficients. The coefficients are found using the 
equality DI = I. 

Example 4. For c = 1 the n-integrals of ^ are 



h — o-xx + bxx ~ + (^xbx ~ b1, 

l[3] 



h = ais] + a-x{bxx - 2a.^x) + 0.% - a^bl- 



Example 5. For c = 2 the n-integrals of ^ are 

Ii = 2a^^ + - 2al + 2ar,b^ - b"^ 



XX 1 

h = a[4] + a^{b[3] - 2a[3]) + a^^{Aa^b^ - 2a\ - h\},\ 

+a'xx{bxx - axx) + b^^a^{a^ - 2b^) + + albl - 2a^6^. 



Example 6. For c = 3 the n-integrals of ^ 



are 



h = flxx + l^xx - (a^)^ + a^bx - ^bl, 

h = a.[6] - 2a[5]a:j.. + bi^la^ + 32a[4]a^ - 30a[4]axbx + lla[4]6^ - A0a[4]aa:x- 
-lla[4]&xx + 146[4]a^ - 15b[4]a^b^ + Yb[i]bl - 106[4]azx - y&[4]&xx+ 



+l9aL + -^bL + 16a[3]6[3] - 36ais]a^^a^ + ISa^s^a^^b^ + SOap]^^ 



-4:5ai3]b^^b^ - 526[3]a^^a^ + 336[3]a^^6^ - 5b[3]b^^a^ - 64a[3]a^+ 
+102a[3]a^5^ _ 2a[3]a^6^ + 13a[3]6^ + 326[3]a^ - 586[3]a^6x' + 386[3]ai.6^- 
-fh]bl + 66aL + fblx - ^^alJ)^^ - ba^^^ + SOalX - ISal^aJ^- 

2'^Xxbx ^^^XxbxX^x ~^ ^'^^XxbxX^xbx "^(^Xxbxxbx '^bxx^xf^X ~l~ SO-ajO^; 

-2Aaxxalbx + 2^axxalbl - Qa^xa^bl + a^xbl - b^^a'^. + ^b^x^lbx - 8bxxalbl+ 
+2hxxaxbl - 2al + Qa% - fa^bl + Salbl - \albl. 

It is remarkable that n— integrals for the differential-difference chain ([7]) given in Examples 
4-6 coincide with ?/— integrals for the corresponding exponential type PDE ([T]) with the same 
Cartan matrix. Recall that in the continuous case these integrals except the last one in Example 
6 were found in [7] (see also [9]). 



4 Integrable cutting of constraints and finite systems 

This section dwells upon alternative approach to construct integrable differential- difference 
chains of exponential type. The approach is based on finding the integrable boundary conditions 
for soliton systems [I9], |21]. In what follows we write functions, that depend on a continuous 
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variable x and two discrete variables n, j as f{x,n,j) = /" and denote its shifts as f{x,n,j + 
k) = fk - Consider the equation 

<-<_^ = e"""'-"?-e""-<^' (26) 

realizing the Backlund transformation for the well known Toda lattice. Here the upper index 
n enumerates the field variables. In equation fl2^ the variables . . . , u~^, u^, v},v?^ . . . are 
considered as dynamical ones. Below we study finite field reduction of the chain. 

Equation (!26|) is the compatibility condition of the following pair of linear equations 



(27) 



n 



which can be written as a discrete version of the Lax equation. Indeed, let U := D„ — u 
V := e"" ~'^'^D~^ + 1, then (l27|l takes the form (f)^ = Ucf), Djcj) = V(f), where Dn, Dj are the 
shift operators acting due to the rule -D„0" = 0"+^, Djcp = Now evidently the compatibility 
condition of the system reads as 

V^ = D,{U)V-VU. 

The latter coincides with the discrete Lax equation. We refer below to fl27|) as a Lax pair. 

Excluding 0""^ from a pair of the equations fl271) we obtain a hyperbolic type linear 
differential-difference equation 



+ ul,<P'l-K + e-''-<^')r = 0. (28) 



The equation fl2B]) is invariant with respect to the transformation defined as n — t- —u, x — ?■ —x, 
n — > —n and j — > —j. Under this transformation the Lax pair (1271) transforms to a Lax pair 



and equation (128|) transforms to the equation 

- C + [-<,. - e"""-""]V^r + = 0. (30) 



(29) 



Thus we have two families fl25]) . fl5Ul) of linear differential-difference equations enumerated 
by n. We would like to know when an equation from one family can be related, by a linear 
transformation, to an equation from the other family. To this end we evaluate the Laplace 
invariants of these equations (see [H]). 

Recall that the Laplace invariants of the equation [T8] : 

y^{3 + 1) + a(j>,.(j) + KjXj + 1) + c{j)y{j) = (31) 

c(j) - a{j)b{j - 1) c(j) - a{j)b{j) - a^{j) 

are given by Ki = and A 2 = • By virtue of these 

au) a(j) 

formulas the invariants Ki^,K2^ and Ki^,K2^ of equations f l28|) . fl30|) are, respectively, 

J^i<t> = e 1 , = e 1 , (d2j 
Ai^ = e -1 ,A2^ = e 1 . (ddj 

It is known that two linear hyperbolic type equations are related to one another by a 
linear transformation only if their corresponding Laplace invariants are equal. Evidently in 
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generic case coincidence of the Laplace invariants generates two constraints on the field variables 
u"'{x,j). Only for some special cases it gives only one constraint. We are interested in such 
special cases. For instance pair of equations Ki^{n,j) = Ki^{n + 1, j + 1), K2^{n,j) = K2^{n + 
1, J + 1) is equivalent to the constraint 



u'^Ti - = - Vj G Z, (34) 

which is interpreted as a cutting off boundary condition for the chain f l26|) . For simplicity we 
put n = 0, so the boundary condition becomes 

uZ\ -u^ = u^ - u\. (35) 

Following [19j we can construct a Lax pair of the corresponding reduction. Under the boundary 
condition fl35l) we have equality of the invariants 



Ki,^(l,j) = i^i,V.(0,j), 

K2,<^(l,j) = i^2,^(0,j) 

and 

7^1,^(1, j) = i^i,0(O, J -1), 
K2,^(l,j) = i^2,<^(0,J-l) 

and we can relate the eigenf unctions 



(36) 



(37) 



(^^ = A{x,3)^\ and 7/;° = S(x,i)0\ (3^ 
Hence, we have from (|271). (|29l). (|38|) 



01 ^ ^ 



(39) 



We study the finite reductions of the chain (!26l) on a finite interval Nj^ < n < Nr. The 
reduction is obtained by imposing the boundary conditions at the left end-point n = 

and respectively at the right end-point n = Nji 

uW=^Nn^^Nn-l_^N^-2_ (41) 

First we concentrate on the left end-point. Due to the reasonings above the eigenf unctions 
should satisfy the following gluing conditions 

(f,^L=A^N,+i and ^^^=50^^+^ (42) 
These conditions allow one to close the Lax equations at the left end-point 

^f^+i = -50^^+1 +M^+V'^"'"'- (43) 
The compatibility conditions of these equations with their counterparts (1^ . (12^ : 
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yield an overdetermined system of equations for A, B: 

A{j + l) = A{j), 



,^L^u^L+^\ AU^ n (45) 

and 

B{i + l) = B{i), 

,iv^+i ~L+n . (46) 



AM + {u^^^ + + e"'^^-"^'^) A{j) = 

Bij + 1) = BU), 

The compatibility of the above equations is guaranteed by the boundary condition (135|) . The 
functions A and B are determined up to a constant and can be chosen so that A - B = —1. So, 
we put 

B = XB and A = -{XB)'^ (47) 

where 5 is a particular solution of (H6l) . A is a constant parameter. 
Put p^{x,j) = ilj^{x,j + 1) for the convenience and find 



'■"~'^-u" in-1 , in J,n _ in+l _ ^,nin 



0n _ g« ^ _ ^n+l _ ^n^n A^L + 2 < 72 < A^^ - 2 (48) 

and 

= e""-""^^"''' + P", p" = + <,,p", for NL + 2<n<NR^2. (49) 

To derive similar equations at the point we use the right end-point constraint fHT]) . for 
which we have 

where 

i(i + l) = i(j), 

Mj) = + + e-f-^-'^^) A{j) = 

and 

S(i + 1) = S(i), 

4(J) = - {u-^ + + e-r^-'^--^) B(,). ^''^ 

The functions A and S are determined up to a constant and can be chosen so that A - B = —1. 
We put A = XAi and B = — (AAi)~\ where Ai is a particular solution of (!50!) . 

And finally, we find in addition to (1481) . (149|) a part of Lax equations corresponding to the 
right end-point 



^^_VH-i ^ e"^ yl0^^«-i+p^^«'\ pi^«-^ = -p^^«-^ + • (52) 

Theorem 4. Chain ^2E) reduced to the finite interval [Nl,Nr\ with the boundary condition 
(2^, is the compatibility condition of the overdetermined system of linear equations 
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Proof. We check the compatibihty conditions of the equations fH5]) - fl52]) for (j)"-, p"-, + 1 < 
n < Nji — 1. It is easy to see that the compatibihty condition for 0^^+^ or for p^^'^^ leads to 

Similarly, the compatibility condition for or for p^'^'^ leads to 

The compatibility condition for 0" or for p" where Nl + 2 < n < — 2 leads to 

u2 - = e"""'-"? - e""-"i^' , NL + 2<n<Nn-2. 
The compatibility condition for A or for B leads to 

And finally, the compatibility condition for A or for B leads to 



In addition to boundary condition fl5Sl) we introduce the degenerate boundary condition 
e~" = 0, for some N and for all integer j. The degenerate boundary condition implies that 
the corresponding eigenfunction is zero, 0^+^ = 0. 

The chain (126!) can be related to the chain 

r--<_^ = e'-""-'^"-^i+''"^\ (53) 

which was considered in the paper [T8|, by the transformation 

m" = r"-i - r". (54) 

The correspondence between boundary conditions is as follows. The condition u^^ = Ui+u\—U2 
leads to the condition 

= rl, (55) 
the condition u^~^^ = —u^^'^ + u^~^ + leads to the condition 

rf = (56) 

and the degenerate conditions e"° = and e~"^^^ = lead respectively to the conditions 

r° = 0, (57) 

r^+i = 0. (58) 
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4.1 Integrable systems corresponding to the Cartan matrices C 



N 



Now we consider reductions of chain (!26l) . We impose condition ( 135|) for n = and the 



degenerate boundary condition e " 



JV+l 



for n = N + 1. The resulting reduction is as follows 



u 



U 



N-1 



U 



N-1 
l,x 



AT— 2 iV— 1 

g„iv-i_„iv 



(59) 



The above reduction admits a Lax pair. Indeed, let us consider vectors > 



p ) and denote 7" = e" "1 . We introduce N x N matrices 



U22 







1 1 














1 





Uu = 





















7^-2 






V 







( 1 












\ 





1 


1' 

















1 




\o 










1 


/ 




( 


-B 


.. 


• M 











.. 


. 




V21 = 







.. 


. 











.. 


• 0^ 









7 



1 

N-1 



( A70 




1 / 



V 



1 

-ul 





1 














\ 







V 

/7° + <x. 

-1 



V22 



\ 













1 


... 





-< 


... 


\ 




«L ... 







-1 u% ... 







-1 


<{3) ) 





where A, B are particular solutions of (H5|) . (H6l) . 

The following lemma gives the Lax pair for system (1591) . 

Lemma 5. The compatibility conditions for the equations (J^j (J^j 



DE 

D-^E 



P 



Uu U12 

U21 U22 



P 



Px 



Vu Vi2 
V2I V22 



P 



(60) 



where U21, V12 are matrices of dimension N with zero entries , E is the unity matrix of dimension 
N, lead to the system (E^j- 



If we impose boundary conditions r ^ 
of chain fl53l) 



l5 



Af+1 



^0 ^0 



l,a; 



'x ' l,x 



r 



r°-r1+2rl 
r-0— r^— rj+r? 



we obtain the following reduction 



r 



N 



l,x 



-2_,.iV-l_,,iV-l+^iV 



„N-1 „JV „]V 



(61) 



which is related to fl^ . This system corresponds to the Cartan matrices C^+i (-Bjv+i). 
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Below we use some manipulations to reduce it to the form fl59p . First rewrite the system as 
follows 



-1 _ ^1 



r ^ = r 



1) 



_ „0 _ „-r"-rV+2rJ 

X l,x 



'x ' l,x — ^ J 



N-1 _ ^N-l _ r'^-2_rN-l_^N-i^M 

I X ' 1 X — ^ ' 

' X ' l,x ^ ' 

fN+l ^ Q 



(62) 



Since r = r} we have — r-^]. = "^^ ^1+^2^ Since r^+-^ = we have r^^^ — r^^^ = 0. 
Thus we can extend the above system to 



J, 1 — ^1 _ prj r|[ r-2+r| 
' a; ' ; 

^0 _ „0 _ -r0-r0+2ri 

'x 'l,a: — ^ 5 

^1 _ ^1 _ „r°-r^-rl+rf 

' X ' l,x — ^ 1 

(63) 

7V_1 _ N-1 ^ .iv-2_^iv-i_^iv-i_^^iv 



N-1 



We subtract the second equation from the first one, the third equation from the second one, 
and so on. Then introducing variables ( 15^ for n = 0, 1, 2, . . . , we arrive at system (!59|) . 



4.2 Integrable systems corresponding to the Cartan matrices A 



N 



We can impose the degenerate boundary conditions at both ends of the chain fl26|) . We put 



and e" 



,N+\ 



that gives the following reduction 



'^x ^l,x 
"^x ~ '^l,x 



U 



U 



N-l 

X 

N 



U 



U 

N 



Af-1 
l,x 



(64) 



The system f l64p admits a Lax pair that can be obtained from the Lax pair f l27p . The 



boundary conditions imply that the corresponding eigenf unctions are zero. 



and 



0. Consider vectors 



We introduce N x N matrices 



R 





1 











... 0\ 




7^ 


1 








... 









1 





... 














1 


V 













7^-^ l) 



( 



and S 





1 





V 











kN+1 



\ 









N-1 



1 



-u 



N 



Lemma 6. The compatibility conditions for the equations 

01 = R(f), (j)x = S(j) 

lead to the system (U^. 



(65) 
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Imposing boundary conditions fl57p and (1581) we get the following chain related with the 
chain (EH) 



' X ' l,x ^ ' 

^2 _ 2 _ T^-r'^-rl+rl 

' X ' l,x ~ ) 



7V_1 7V-1 _ .iv-2_^iv-i_^iv-i_^^iv 
' X ' \,x — ^ ■ 

l,x 



(66) 



The above system corresponds to the Cartan matrices A^. 
The system ( 166|) can be reduced to the system 



1 1 2 

II-'- n-"- . 11-^ 



- vl^^ = -e" - e"! + e'^i, 
t;2-v2'^ = -e'^' -e''? +6^^?, 



IV_1 AT-l t,JV-2 iV-1 N-1 N 

'^ x ~ '^ Ix = 6 ~ ^ ~ ^ + e > 



(67) 



■u" — vi\. = e — e — e 1 



called two-dimensional Toda molecule equation, found in paper [13], by the transformation 

e"" = r: - rl,. (68) 



The correspondence between boundary conditions is as follows. The conditions r° = 0, r 



N+l 



lead respectively to the conditions e^" = 0, e""^^^ = 0. 

(2) 

4.3 Integrable systems corresponding to the matrices Djy 

We can impose boundary conditions ( 1M|) at both ends of the chain (12B]1 . That is for n = we 
have (1351) and for n = N we have 

This leads to the following reduction 



(70) 



Imposing boundary conditions (155|) and (l56ll we get the following chain: 



j.0_j.0_,_2j.1 



/>^0 /v>0 — 

' x l,x ~ 



' X 1 X — ^ ' 

'x — ^ ) 



(71) 



which is related to (170|) . This system corresponds to the Cartan matrices -D^^^ 



17 



Consider vectors 
matrices 



. . . , (j)^-^),p = (p\ p^ . . . , p^-i). We introduce A^- 1 x iV - 1 



Uii 



/ 1 

7^ 1 

. 

V . 







/ 





U21 



7^-3 







1 

N-2 



( A70 




1/ 



V 







-Mi 



-ul 



V .. 

1 



^7^-1 y 



1 









V 







V2\ 



/ -5 



\ 





... 

. 0\ 

. 

. 

• 0/ 



—u 



N-2 



V22 



( 1 7' 

1 

... 

Vo ... 

^ 




r 





7 



. 

• oy 

\ 


N-2 








1 / 








... 


\ 







... 





V\2 = 










u 


... 


B ) 



( 7i + ul 
-1 




V 
























-1 


"^1,1: 










-1 u^' 



(J) / 



where A, B, A, B are particular solutions of (SSl), 1^, ([50]), (JST 
The following lemma gives the Lax pair for system (ITOl) . 



Lemma 7. The compatibility conditions for the equations [4^ ), (J^j ^5^) . (5l\) . 



DE 

D-^E 



P 



Un U12 
U21 U22 



P 



Pa 



VU V12 
V2I V22 



(72) 



where E is the unity matrix of dimension N — 1, lead to the system ((70| j 
Example 7. Put N = 3. That is we have the chain 



u 



n = 0, 1,2,3, 



(73) 



with boundary conditions u ^ ~v!l = u\—u\ and u^i—u^ = — u\. So we have the system 





= u\- 


ul 


< 






< 



















'4 — e"^" 



(74) 



2_„,3 ,,1 _ 

1 _ e"-i 



Let us give the Lax pair for the system in an explicit form: the compatibility conditions 
for the equations (2^, JT^ ; ( f50l) . ( f37]) and 



DE 

D-^E 



Mil U12 
U2I U22 



Px 



Vll V12 
V2I V22 



(75) 
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where (p = {(j)^, 0^)^, p = (p^, P^)'^, 



10\ f A-f° 0\ f \ [ I 1^ 

7^ 1 j ' = I J ' = U ' ^ V 1 



-ul 1 





f 




) , ^12 = 1 


V 


B j 



// we impose boundary conditions = tJ, = we obtain the following reduction of 
chain ( f3^) 

ri-rj, = e^"-^'-^i+'-^ (76) 



^0 _ ^0 _ -r°-r°+2r\ 



^2 _ 2 _ 2ri-r2-r2 



(2) 

The above system corresponds to the Cartan matrix of . Below we use some manipulations 
to reduce it to the form IjiTW. First rewrite the system as follows 



„0 ^0 



1 

.0 «0_|_«1 



„2 ^2 



Since r~^ = r{ we /iai'e r~^ — rj~^ = e'^°~'^^~'^^'^^^ . Since r^ = r'^^ we have — rf^ 



g?-_i r +r ^ j'fmg yj^ gan cxtcnd the above system to 



-1 ^-1 



I 



rl-rl-rl+r^ 



^0 _ ^0 _ r-l-rO-rO+rl 



r^-ri,. = e^"-^-'-J+'^?, (78) 

'x ' l,x — ^ 

^1 „1 



We subtract the second equation from the first one, the third equation from the second one, and 
so on. Then introducing variables ( [5^ for tt, = 0, 1, 2, 3 we arrive at system ( f^] ). 



5 Conclusions 

Integrable differential-difference exponential type systems are presented corresponding to any 
simple or affine Lie algebra. In fact they are discrete analogues of the well-known Toda field 
equations. In the literature such kind systems are called lattice Toda field equations (see, for 
instance, [TB], [IZ]). Our systems are different from those found earlier in |17j except the case 
Ajq. The system corresponding the algebra Ajsi called two-dimensional Toda molecule equation 
has been found in [T3] . 

We suggested a formal procedure which assigns to any Cartan matrix a system of differential- 
difference equations. Actually this correspondence has an important algebraic interpretation. 
Recently it was observed that any differential-difference system of the form 

= /*(u, u^., ui), j = 1, A;, u = (m\ u''), u = u(n, x), ui = u(n + 1, x) 



19 



admits characteristic Lie algebras in x and n directions. It is shown in the article that the 
characteristic x-algebras of the systems ([3]) are isomorphic to the Borel subalgebras of the 
simple Lie algebras related to the corresponding Cartan matrices. 
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